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q \ We study two different decoherence modes for entangled qubits by considering a 

Liouville - von Neumann master equation. Mode A is determined by projection 
^ | operators onto the eigenstates of the Hamiltonian and mode B by projectors onto 

Q \ rotated states. We present solutions for general and for Bell diagonal states and 

çf) ' calculate for the later the mixedness and the amount of entanglement given by the 

concurrence. 

We propose a realization of the decoherence modes within neutron interferometry 
by applying fluctuating magnètic fields. An experimental test of the Kraus operator 
decomposition describing the evolution of the system for each mode is presented. 
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3 ! I. INTRODUCTION 

cr 

Closed quantum systems are idealizations which do not exist in a real physical world. 
Actually, one always has to deal with open quantum systems which arise due to an 
interaction of the system under consideration with an external environment (e.g. reservoir 



heat bath) [JQQ. The system - environment interaction causes a phenomenon known 
as decoherence: quantum correlations and interferences are destroyed in course of time; 
the system shows more and more classical behavior. The theory of decoherence is one 
candidate to solve the question why our world looks so classical 0, ■ 

The total Hamiltonian of system and environment generates a unitary time evolution 
U (i) and is of the form Hse = H®t E + t^H E -\-Hj, where H, H E and Hj are, respectively, 
the system, environment and interaction Hamiltonians. The evolution of the system, 
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represented by the density matrix p(t), or the reduced dynamics is obtained by tracing 
over the environmental degrees of freedom p(t) = Tte ps+EÍf) = TíE{U{t)ps+EÍ^)U\t)) 
and thus inheriting a nonunitary evolution for the system in contrast to closed systems. 

In most of the cases we do not have access or information about the dynamics of the 
environment. Therefore we have to describe the evolution of the system by an effective 
dynamics: the Liouville - von Neumann master equation. Thereby it is not so important 
to know the exact Hamiltonian and the nature of the environment but only its effects on 
the system. Our strategy in this paper is to propose several effective models which do not 
care about the exact nature of decoherence but provide scenarios how decoherence can 
affect a system. 

Under several assumptions [ïj], such as Markovian semigroup approach, complete posi- 
tivity, initial decoupling of system and environment, and weak coupling, the dynamics of 
the system can be described by a Liouville - von Neumann master equation 

^p(t)=-i[H(t),p(t)]-D[p(t)]. (1) 

Lindblad and Gorini-Kossakowski-Sudarshan 0, 0] derived the most general structure of 
the dissipator 

D[p(t)] = ^(A[4p(í) + p(t)A{A k - 2A kP (t)A{) , (2) 

k 

where A k represents a so-called Lindblad generator. The sum is taken over an arbitrary 
number of components but maximally up to n 2 — 1, where n denotes the dimension of the 
system. For simplicity we choose the generators to be projectors such that A k = y/\k~Pk 
with Pj 2 = Pk (see Ref. 0) which gives for the dissipator 

D \P\ = \J2 Xk ( Pk P + P Pk ~ 2P kP p k) • (3) 

k 

The paper is organized as follows. In the next section we introduce and discuss two 
possible decoherence scenarios for a two qubit system by choosing different projection 
operators Pk- In Sect lHIl the two decoherence modes are discussed for the special case of 
Bell diagonal states. In Sect lIVI we propose a realization of the decoherence modes within 
neutron interferometry via random magnètic fields which represent the environment. In 
Sect]V]we present the Kraus operator decomposition. The action of this decomposition is 
mathematically equivalent to the Lindblad form of the Liouville - von Neumann equation. 
We can test this equivalence by a simple experiment with single neutrons. 



II. DECOHERENCE MODES IN A TWO QUBIT SYSTEM 

Let us consider a two qubit system with Hilbertspace H = Tv- 1 ' ®üS 2 ' = C 2 ®C 2 where 
{| e fc)}fe=i,...,4 denotes an eigenbasis defined by H\ek) = Ek\ek), with H = H^0t + t®H^ 
the Hamiltonian of the undisturbed system. A general state p of the system can be 



expressed in the eigenbasis p = J2k j Pkj\ e k)(zj\, where (pkj) denotes the 4x4 coefíicient 
matrix. 

We consider Lindblad generators Pk that project onto one-dimensional subspaces and 
fulfil Pk = 1, furthermore we assume that only one dissipation parameter À parameter- 
izes the strength of the interaction and therefore of the decoherence. Then the dissipator, 
Eq.flHJ), can be written as 

4 

D[p} = X{p-J2PkpPk) ■ (4) 

fe=i 

In the following sections we solve the Liouville - von Neumann equation with the 
dissipator (jlj) by assuming different projection operators Pk, what we call decoherence 
modes. 



A. Mode A 



The first mode describes the simplest possible case. The Lindblad generators are chosen 
to be projectors Pk = |efc)(e&| onto the eigenbasis of the Hamiltonian (see Refs.jaBl)- In 
this mode of decoherence the time evolution (0) for the coefficient matrix is given by 

Pkj = (-i{E k - Ej) - X A ) pkj for k ^ j 

Pkk = , 

which can be easily solved 

Pkj(t) = e-^-^e^pkM for k ± j 
Pkk{t) = pfcfc(O) . 

The decoherence affects only the off-diagonal elements and leaves the diagonal elements 
untouched. 



B. Mode B 

For the second mode the Lindblad generators are chosen to be projectors Pk = \ek)(èk\ 
onto the following states 

v 1 K 1 

|ei,3) = -^=(|ei)±|e 3 )) , |e 2 , 4 ) = ^(|e 2 ) ± |e 4 )) , (7) 

where the upper (lower) sign corresponds to the first (second) index. 

The time evolution of the coefficient matrix can be separated into 3 types of differ- 
ential equations. Type I is vàlid for the components p í2 , pu, P23, P34 and has the structure 



P12 = (-i{Eí - E 2 ) - \ B )pi2 , 



(8) 



in analogy to mode A. Type II holds for the diagonal components and reveals pairwise 
coupled differential equations for pn - p 33 and p 2 2 - Pu of the form 

Xb , Ab . Xb Xb 

Pu = — 2~Pii + "yPss , P33 = -yPii - -ypss • (9) 

Type III also gives pairwise coupled differential equations 

PÏ3 = {~i{E l - E 3 ) - -y)pi3 + ^Psi , P31 = ^Pia + (i{Ei - E 3 ) - y-)p3i , (10) 

vàlid for the components P13 - P31 and P24 - P42. The solutions for the several types of 
differential equations are: 
for type I, Eq.flBJ), 

Pi2(í)=e-^- £2 V A *Vi2(0), (11) 

for type II, Eq.flÜJ), 



Pn(í) = ~(1 + e- Aflí )pn(0) + i(l - e- Aflí ) P 33(0) , 
M*) = \{l ~ e- ABÍ )pn(0) + 1(1 + e^)p 33 (0) , 
for type III, Eq. ífTüjl . 

p 13 (t) = e-^((coshf - fc^sinhf )p 13 (0) + ^ sinh £ p 31 (0) 
p 31 (t) = e-^ ((cosh f + smh f )p 31 (0) + ^ smh f p 13 (0) 



(12) 



(13) 



where p = ^/A| - A{E 1 - £ 3 ) 2 . 

Decoherence mode B affects not only the off-diagonal elements of the density matrix 
but also the diagonal ones. 

In the case we consider in this paper - decoherence modes in neutron interferometry, 
Sect JIVI - the system is given by the free neutron passing through the interferometer, 
whereas the magnètic fields placed in represent the external environment. Thus there is 
no splitting in the energies, e.g., E\ = E 3 , E 2 = £4 and p = Xb, so that type III is equal 
to type II. 

Remark. It is worth noting here that the choice of projection states, Eq. (|7|) . corresponds 
to a rotation of states in one subspace. Suppose we split the eigenstates of the undisturbed 
Hamiltonian in eigenstates of the subspace Hamiltonians {|ai), \a>2)} and {\bi), I&2)} m the 
following way 

|ei, 3 ) = |oi,2)|&i) , |e 2 ,4> = \ai,2)\b 2 ) • (14) 



Now consider a rotation of the first subbasis, {|+) = -^(|ai) + |a2)), |— ) = ^(\ a i) — \ a 2))}, 
the second subbasis is left untouched. The basis of the total Hilbertspace changes 



|ei) = |+>IM = ^(|ei> + |e 3 », 
\e 3 ) = Hlb,} = ^=(\ ei ) - \e 3 )) , 



\e 2 ) = \+)\b 2 ) = ^=(\e 2 ) + \e 4 )) 1 
\eí) = \-)\b 2 ) = ^(N) - l e 4)) , 



(15) 



which corresponds exactly to the states used in Eq.®. Therefore decoherence 
mode B can be denoted "R <g> E " to indicate the rotation of the first subspace and the 
untouched eigenbasis in the second subspace whereas mode A can be labelled by u E<g>E" . 

In the case of photons (see, e.g., [ï(j) the eigenbasis E corresponds to horizontal \H) 
and vertical \V) polarization whereas the rotated basis R represents polarization states 
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0) 



wc 



|+45°) and |— 45°). In the case of neutral kaons (for an overview see, e.g. 
can identify the eigenbasis E with the short- and long-lived states \K$) and \Kl) and the 
rotated basis R with \K°) and \K°). 



III. INITIAL CONDITIONS BELL DIAGONAL STATES 

We want to illustrate the above discussed decoherence modes by choosing a certain 
class of states as initial conditions - the so-called Bell diagonal states p = i) i\ 

with Yli v % — 1) which are diagonal in the Bell basis 



|v]/ li2 > = — (| ei > ± |e 4 » , 
In the Standard basis they are expressed by 



l*3,4> = -^(|e 2 )±|e 3 }) 



(16) 



fvi + v 2 



p 











v x - v 2 \ 



z/ 3 + z/ 4 z/ 3 — z/ 4 
u 3 — v± z/ 3 + z/ 4 
\v x -v 2 + v 2 ) 



/Ei Ai\ 
S 2 A 2 
A 2 S 2 
\Ai 



(17) 



using the notation Ej = v\ + v 2) E 2 = u 3 + z/ 4 , Aj = v\ — u 2 , A 2 = u 3 — v±. 

States are characterized by two quantities: mixing and entanglement. The mixedness, 
defined as 5 = Trp 2 , ranges between 1 (pure states) and \ (maximally mixed states) 
and is given by 5 = v\ + v\ + v\ + v\ for Bell diagonal states. The concurrence C 
[Üï- JJ, 3| is a suitable quantity that measures the entanglement contained in a state p. 
It is defined by C(p) = max{0, p\ — p 2 — p 3 — /i 4 }, where pi are the square roots of the 
eigenvalues in decreasing order of the matrix R = p (a y ® a y ) p* (a y <8> <? y ) and p* denotes 
complex conjugation in the Standard basis. The concurrence varies between 1 (maximally 
entangled states) and (separable states) and is given by C = max{0, 2max{z/j} — 1} for 
Bell diagonal states, depending on which weight z/j is the largest. A Bell diagonal state 
can only be entangled (C > 0) if the largest eigenvalue fulfills v^>\. 



The special case of a pure and maximally entangled Bell state, e.g., the Bell singlet 
state 1^/4), where U4 — I and v\ = v 2 = ^3 = or £1 = Ai = and £ 2 = — A 2 = 1, yields 
5 = 1 and C = 1. 



Mode A 



The initial Bell diagonal state, Eq.ífTT)). evolves in time according to mode A, Eq.( 
into the state 





/ 


£1 










1 







£2 


e~ AA *A 2 





2 







e~ XAt A 2 


£ 2 

















Ei / 



;is) 



For the mixedness of the state we get 6 = |(£? + £ 2 + e 2A a*(A 2 + A 2 ,)) and for the 
concurrence we find C(p) = max{0, 2 max{/ij} — l}, where /íi i2 = §(£1 ± e~ AAÍ A!) and 
/ i3,4 = |(S 2 ±e^ í A 2 ). 

Choosing the Bell singlet state ^4) the mixedness ó = |(l + e~ 2AAÍ ) ranges from a pure 
state (5 = 1) to a mixed but not maximally mixed state (5 t ^°°> |). The concurrence 



C(p) = e AÍ decreases exponentially from a maximally entangled state (C = 1) to an 
asymptotically separable state (C f ^°°> 0). The behavior of 5 and C is plotted in Fig^ 



B. Mode B 



The second mode, Eqs.ffTT |) -(fTïï j) . generates the density matrix 
/l 



e^ ABÍ A 



Pit) = Z 
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(19) 



with the notation A = £1 

We obtain for the mixedness 8 = \ (l + e~ 2Así (2A? + 2A 2 . + (£1 + £ 2 ) 2 )) and for the 
entanglement C(p) = max{0, 2 max{/ij} — l}, where /íi j2 = 4(1 + e~ Así (A ± 2A 2 )) and 
/i 3 ,4 = K 1 - e " Aflí ( A T2A 1 )). 

The special case of ^4) yields the following results. The mixedness 5 = 4(1 + 3e~ 2Aflí ) 



1) to a maximally mixed state (S 



t—*oo 



The concurrence 



varies from a pure state (5 
C(p) = maxjO, 7j(3e~ Así — 1)} decreases exponentially and the initially maximally entan- 
gled state (C = 1) reaches the border of separability (C = 0) at finite time t — where 
the mixing has the value of S — | (see 
with respect to Xt is shown. 



also Ref. [3]). In Fig|T]the dependence of 5 and C 
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FIG. 1: Graphical comparison of the mixedness 5 and the concurrence C in dependence of Xt 
for mode A and B for the Bell singlet state 1^4). The upper curves correspond to mode A and 
the lower ones to mode B. 



IV. REALIZATION OF DECOHERENCE MODES FOR NEUTRON STATES 



The different decoherence modes presented in Sectjn]and discussed for Bell diagonal 
states in Sect IIIIl can be tested within neutron interferometry 17]. A neutron is entangled 
18 , IÏ9I EoT ] between the internal degree of freedom - spin - and the external degree of 
freedom - path - which is described by the bipartite Hilbertspace H = 7i spin ®H pa th- Let 
us consider the antisymmetric Bell state which is experimentally feasible 



l^exp) = |*4> = ^=(|ft> ® [H> 



® |I)) = -j={\e 2 ) - Ics» , 



(20) 



where \ f) and represent ±z polarized spin states whereas |I) and |II) denote the paths 
in the interferometer. 

Calculating the density matrices for both modes A and B we find 



p A {t) = - 



/o 0\ 

1 - e - XAt 

-e~ XAt 1 

\0 0/ 



(21) 



p B (t) 
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n 
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-X B t 











1 + e _As * - 
-2e~ XBt 1 






-2e~ XBt 
+ e" As * 










-,— A_BÍ 



(22) 



The off-diagonal elements fade away exponentially for both modes. For mode B the 
diagonal elements are distributed in the whole 4-dimensional space so that at t — > 00 the 
density matrix approaches the normed unity, i.e., the totally mixed state. 

Within neutron interferometry all matrix elements can be determined experimentally 
[2~l·l | via the procedure of quantum state tomography ET 



A. Decoherence via random magnètic fields 



For the implementation of decoherence we use randomly fluctuating magnètic fields 
which act on an ensemble of neutrons produced in the specific state p. 

The action of a magnètic field B = Bn in the direction n on a neutron state is described 
by the unitary operator U(a) = e*^™' ", where a = 2psBt = uj^t denotes the rotation angle 
and pb-i &l the Bohr magneton and Larmor frequency, respectively. 

The neutron beam passes a fluctuating magnètic field in such a way that each neutron 
which is part of the quantum mechanical ensemble described by p feels separately a 
different but constant magnètic field. This corresponds to applying a unitary operator 
U(a) with constant rotation angle a onto the density matrix p. For the whole ensemble 
we have to take the integral over all possible rotation angles a 

p — >pf = J U(a)pU\a\ P{a)da , (23) 

p(a) 

where P(ot) denotes a distribution function. In our case the distribution function is a 

Gaussian Pia) = r— e ^ with Standard deviation a. Although each transformation 
separately is unitary due to the integration we end up with a nonunitary evolution. 



B. Mode A 

For an incoming polarized neutron the state l^exp) is prepared after passing the beam 
splitter and spin flipper. This initial state is subjected to the fluctuating magnètic fields 
oriented along the z-axis in each path of the interferometer, see Fig|21 The rotations 
U(a) and U(/3) caused by the fields are independent but their distributions have the same 
deviation o. 

The action of the two magnètic fields can be described by a "conditioned operation". 
Depending on the state of the spatial degree of freedom either operation U(a) or U((3) is 
applied to the spin state 



|^ S pin) ® |I> — > U(a)\ifj spin ) (8) |I) 

IVspin)® |H) — »• U{P)\ ^ S pi n }®|II) 



(24) 



For a single neutron the application of the conditioned operation on the initial state |^/ e xp)j 
Ea. ípüj) . gives 



p(a,{3) = - 
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(25) 



Spin flipper 




B z ll \(3) 



Phase shifter x 



FIG. 2: Experimental setup for the realization of mode A. The spin flipper is inserted to achieve 
the Bell singlet state. The magnètic fields Bz (ot) and Bf l \(3) produce independent rotations 
U(a) and U((3), respectively. With the phase shifter \ an d the spin rotator £ the final state is 
analyzed. 



which after integration over a and f3 turns into 



p = / p{a,(3)P{a)P{(3)dad(3 
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By comparison of Eq. ÍI21I) and Eq. (|2fi|) we immediately see that 



(26) 



(27) 



the decoherence parameter is directly related to the deviation a of the fluctuating 
magnètic fields. Note, that for only one magnètic field located in one of the paths fluctu- 

2 

ating with deviation a the above relation is given by X^t = and for one field acting 

2 

on both paths we have Xaí = \. 



C. Mode B 

For mode B we prepare the same state l^exp) but use different fluctuating magnètic 
fields, shown in Fig|ïï| The different unitary operations caused by the magnètic fields are 
assumed to act independently but the Gaussian distributions have the same deviation a. 

In order to implement experimentally the rotated projectors (jTJ) we need a magnètic 
field in x-direction B x . However, to achieve the same damping in the off-diagonal elements 
as in the diagonal elements we have to insert an additional magnètic field in z-direction 
B z which influences only the off-diagonal elements. It reflects somehow the effect of the 
Kraus operators which act in x- and z-direction (see operator M3 in Eq. (jH4j) ). 
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FIG. 3: Experimental setup for the realization of mode B. The magnètic fields B% (7), -Bi llj (ó) 



Si (a) and B^ L> {(3) generate independent unitary rotations Í7(7), í/(a) and U((3), re- 

spectively. The order of the magnètic fields in each path does not matter in this context. 



>(H) 



The neutron after the conditioned operation is described by the density matrix 
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For the ensemble state we get after the integrations over the angles a, (3, 7, 6 with 
Gaussian weights 



/l 
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(29) 



and by comparing Eqs. (j22j) and ()29|) we find the relation 

■2 



X B t = — 
2 



(30) 



between the decoherence parameter A# and the deviation a of the Gaussian distribution. 



Experimental test. Our decoherence modes A and B can be tested experiment ally in 
the following way. The incoming polarized neutron, prepared after the beam splitter and 
spin flipper in a Bell singlet state, is subjected to magnètic fields with a certain Gaussian 



variation a in both paths of the interferometer (see FigJ2] for mode A and FigO for mode 
B). Due to relation (|27|) for mode A and (jrU))) for mode B the value of the decoherence 
parameter À - the "dephasing" due to the variation of the magnètic fields - is adjusted. 
The time t corresponds to the duration the neutron remains in the interferometer, more 
precisely, within the magnètic fields and remains constant. The individual density matrix 
elements are measured experiment ally via quantum state tomography and have to be 
compared with the corresponding theoretical expressions (J2~ï|) and í)22|) . By varying a, 
which means varying À, one can nicely examine the specific exponential decrease of the 
decoherence modes A and B. 



V. CONNECTION TO THE KRAUS OPERATOR DECOMPOSITION 

In the following section we present a connection to the Kraus operator decomposition. 
According to the theory of decoherence 0,1^01 the non-unitary evolution of the system can 
be described by Kraus operators. We want to demonstrate that by a simple experiment 
within neutron interferometry. We can check whether the theoretically predicted Kraus 
operators correspond to the implemented decoherence modes discussed in Sect IIVI 



A. Kraus operator decomposition 

The completely positive time evolution generated by the Liouville - von Neumann 
master equation (JTJ together with the Lindblad form of the dissipator (J2J) can also be 
represented by a dynamical map expressed in the Kraus operator decomposition 0, 

p(0)^p(t) = J2MkP(0)Ml, (31) 

fe 

where the Kraus operators M k fulfil M^Mk = 1. The first approach corresponds to a 
continuous time-dependence of the state whereas the second one treats decoherence via 
discrete state changes. Both views are equivalent and a correspondence between Lindblad 
generators A k and Kraus operators M k exists for small 5t (see e.g. Ref . ) 



M = l-(iff + i$>í4b)« (32) 
M k = VótA k . 

Clearly, the Kraus operators are not uniquely determined by Eq. íKUj) and allow for a 
unitary transformation. 



B. Mode A 



For the Hilbert space we are using, H spm <8> 7Y pat h, where entanglement occurs between 
spin and spatial degrees of freedom, mode A represents a kind of phase flip channel 



which destroys the coherence of the system. In this case the Kraus operators are given 
by 



M = Wl - ^ V ® l p Mi = J j l s ® 

V Vi (33) 

where w = Xt is the probability for the occurring decoherence. It leads for small St to the 
state Eq.fjüJ), which allows for a general initial condition. 
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FIG. 4: Experimental setup for the realization of the Kraus operator a s z <Jz for mode A. 

Experimentally the Kraus operators can be implemented in the interferometer in the 
following way Incoming polarized neutrons are prepared as Bell singlet states and feel 
the effect caused by the Kraus operators, as shown in Fig0] for the operator M3. The 
identity operators l s and V clearly do not change the spin and spatial degree of freedom. 
The operator a s z when acting on the spin state |JJ.) induces a phase shift of tt. This phase 
shift difference between spin up and spin down can be implemented by a magnètic field in 
z-direction B z (modulo an overall phase shift). The operator a p z on the spatial subspace 
is realized by a phase shifter in the path |II) which induces a fixed phase shift of n. 

The states produced by the four Kraus operators are measured tomographically and 
the weighted sum according to (jïïïïj) represents the state p, Eq.(jHJ), of mode A. 



Mode B 



Mode B is a combination of a bit flip channel and a phase flip channel |2|, |24|. The 
corresponding Kraus operators are 



y (34) 
M 2 J|<®F Mz=Jjal®al, 



and create for small St the state given by Eqs. (fT2*|) and (|TÏÏ|) (u; = At) allowing for 

general initial conditions. 




Phase shifter (\) 



Br Phase shifter (ix 



FIG. 5: Experimental setup for the realization of the Kraus operator a x ® a z for mode B. 



The Kraus operator for the spatial part a p z is the same as for mode A, inducing a phase 
shift of 7T. The difference to mode A lies in the a x operator for the spin part. It can be 
realized by two magnètic fields B x in both arms pointing in the x-direction, which cause 
a spin flip, see FigJ3J 

Again, the weighted sum of the measured states produced by the Kraus operators 
according to (jMjl leads to the state of mode B. 



VI. SUMMARY AND CONCLUSION 

We have considered the Liouville - von Neumann equation where decoherence is im- 
plemented by the dissipator in Lindblad form. We study two kinds of decoherence modes 
where the Lindblad generators are given by different projection operators: decoherence in 
the eigenbasis of the Hamiltonian, mode A, and decoherence in a rotated basis, mode B. 
The two modes are analyzed in detail for Bell diagonal states, where it turns out that in 
mode B the state gets more mixed and the entanglement decreases faster than in mode A. 
The Bell singlet state \^ À) gets separable at finite Xt = ln 3 in mode B whereas in mode 
A the state remains still entangled at that point by an amount of 33%. 

The realization of the proposed decoherence modes uses the bipartite Hilbert space 
construction of neutron interferometry where entanglement for single neutrons occurs 
between an infernal (spin) and an external (path) degree of freedom. 

We create decoherence via magnètic fields in the interferometer and find that the 
decoherence parameter À is determined by the deviation a of the fluctuating fields, Eqs. 
(12 7|) and (pïïl|) . This allows an experimental control of the implemented decoherence in 
each mode. The strength of decoherence does not depend on the actual rotation parameter 
a of the magnètic field but only on the width of the Gaussian distribution. 



Measuring experiment ally the matrix elements of a state via state tomography and 
varying a we examine the time evolution of the state according to mode A and mode B, 
Eqs. (pHj) and p2j) . 

In addition we can test experiment ally the validity of the Kraus operator decomposition 
which alternatively describes the completely positive time evolution. The Kraus operators 
are constructed for each mode and realized within neutron interferometry. 
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